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ABSTRACT 

Recently published Spitzer Space Telescope observations of the classical 
Cepheid archetype 5 Cephei revealed an extended dusty nebula surrounding this 
star and its hot companion HD 213307. At far infrared wavelengths, the emission 
resembles a bow shock aligned with the direction of space motion of the star, in- 
dicating that 5 Cephei is undergoing mass-loss through a stellar wind. Here we 
report Hi 21-cm line observations with the Very Large Array (VLA) to search 
for neutral atomic hydrogen associated with this wind. Our VLA data reveal a 
spatially extended H I nebula (~ 13' or 1 pc across) surrounding the position of 
S Cephei. The nebula has a head-tail morphology, consistent with circumstel- 
lar ejecta shaped by the interaction between a stellar wind and the interstellar 
medium (ISM). We directly measure a mass of circumstellar atomic hydrogen 
Mhi ~ 0.07 M Q , although the total Hi mass may be larger, depending on the 
fraction of circumstellar material that is hidden by Galactic contamination within 
our band or that is present on angular scales too large to be detected by the VLA. 
It appears that the bulk of the circumstellar gas has originated directly from the 
star, although it may be augmented by material swept from the surrounding 
ISM. The Hi data are consistent with a stellar wind with an outflow velocity 
V = 35.6 ± 1.2 km s" 1 and a mass-loss rate of M « (1.0 ± 0.8) x 10~ 6 M yr _1 . 
We have computed theoretical evolutionary tracks that include mass loss across 
the instability strip and show that a mass-loss rate of this magnitude, sustained 
over the preceding Cepheid lifetime of 5 Cephei, could be sufficient to resolve 
a significant fraction of the discrepancy between the pulsation and evolutionary 
masses for this star. 
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1. Introduction 

Owing to the tight coupling between the pulsation properties of Cepheid variables and 
their fundamental stellar parameters, these stars have long provided crucial tests for stellar 
evolution models (Gautschy & Saio 1996) and play a key role in the calibration of the 
extragalactic distance scale (Feast & Walker 1987; Freedman et al. 2001). Nonetheless, for 
more than 40 years, a puzzle has persisted concerning the lingering discrepancy between 
stellar masses derived for Cepheids using different methods (e.g., Cox 1980). In spite of 
improved radiative opacity calculations (Iglesias et al. 1990; Seaton et al. 1994), advances 
in evolutionary modeling (Bono et al. 2002), and the consideration of metallicity effects 
(Keller & Wood 2006), masses derived from stellar evolutionary models are found to be 
systematically ~10-15% higher than those derived from stellar pulsation models (Caputo et 
al. 2005) or orbital dynamics, when available (Evans 2009). 

As described by Neilson et al. (2011), one of the most promising solutions to the so- 
called "Cepheid mass discrepancy" is likely to be convective core overshooting (Chiosi et al. 
1992) coupled with mass-loss during the Cepheid phase of evolution (see also Bono et al. 
2006). Mass- loss from Cepheids has long been postulated, not only as a means of resolving 
the Cepheid mass discrepancy, but as a natural consequence of stellar evolution and pulsation 
models (e.g., Iben 1974; Willson & Bowen 1984; Neilson & Lester 2008). The observational 
confirmation of mass loss from Cepheids would impact not only our understanding of Cepheid 
evolution, but would also have implications for the use of Cepheids as distance indicators. 
For example, mass-loss can non-negligibly affect the structure and scatter in the infrared 
period-luminosity relation (Neilson et al. 2009), and the presence of circumstellar material 
can impact Cepheid distance determinations made using the interferometric Baade-Wesselink 
method (Merand et al. 2007). 

A number of previous authors have attempted to identify direct evidence of past or 
ongoing mass-loss from Cepheids using observations ranging from ultraviolet to radio wave- 
lengths (e.g., Deasy & Butler 1986; Deasy 1988; Welch & Duric 1988; Bohm-Vitense & Love 
1994; Merand et al. 2007; Kervella et al. 2006, 2009; Neilson et al. 2009). However, the 
results have been inconclusive, as the mass-loss rates (or upper limits) derived from these 
studies (M 10~ 12 to 10~ 5 M yr _1 ) span many orders of magnitude and are typically 
rather uncertain owing to their reliance on trace atomic and molecular species (which may 
not be optimal tracers of past and ongoing mass- loss), their inability to sample very extended 
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spatial scales, and finally, the strong dependencies of derived mass-loss rates on a variety of 
underlying assumptions (e.g., wind structure, gas-to-dust ratios, ionization fractions). 

With the goal of obtaining new empirical constraints on Cepheid mass-loss, Marengo 
et al. (2010a) recently used the Spitzer Space Telescope to survey a sample of 29 nearby 
Cepheids. While IR excesses were not directly detected, ruling out the presence of a large 
amount of warm (~500 K) dust in close proximity to the stars, the Spitzer images revealed 
extended emission around several targets (Marengo et al. 2010a; Barmby et al. 2011). 
Perhaps the most intriguing was the discovery of a prominent nebula around the archetype 
Cepheid variable, 5 Cephei (Marengo et al. 2010b). 

Some basic stellar properties of 5 Cephei are summarized in Table 1. 5 Cephei is part 
of a wide binary (Benedict et al. 2002), with a hot companion HD 213307 (T e g =8800 K; 
Cenarro et al. 2007) separated by a projected distance of 40". The IR emission surrounding 
S Cephei (and HD 213307) is visible in multiple Spitzer bands and shows a roughly parabolic 
structure with an extent of ~ 5' (~ 2 x 10 4 AU). The symmetry axis of the parabola is aligned 
with the direction of the star's motion through the interstellar medium (ISM; Marengo et 
al. 2010b), where V space « 10.3 km s" 1 and P.A.=58°.3E The IR emission therefore appears 
to trace a bow shock structure, as arises when a moving, mass-losing star interacts with the 
local ISM (e.g., Wilkin 1996). Bow shocks have been seen previously in FIR images of a 
number of mass-losing stars, including the supergiants a Cam (Van Buren & McCray 1988) 
and a Ori (Noriega-Crespo et al. 1997; Ueta et al. 2008) and the asymptotic giant branch 
(AGB) stars R Hya (Ueta et al. 2006; Wareing et al. 2006), Mira (Ueta 2008), R Cas (Ueta 
et al. 2010), IRC+10216 (Ladjal et al. 2010), TX Psc, and X Her (Jorissen et al. 2011). The 
emission mechanism responsible for the FIR emission from these circumstellar bow shocks is 
uncertain, but is likely to be mainly thermal emission from dust, with possible contributions 
from low-excitation atomic emission lines (Ueta et al. 2008; Marengo et al. 2010b). 

While the companion of 5 Cephei, HD 213307, cannot strictly be excluded as the source 
of the IR emitting material seen by Spitzer, this appears quite unlikely based on the B7- 
B8 III-V spectral type of this star (see Marengo et al. 2010b for discussion). Thus the Spitzer 
observations of Marengo et al. (2010b) provide strong evidence that 5 Cephei is undergoing 
mass loss through a stellar wind. To further characterize the nature of this wind and its 
interaction with the ISM, we have now used the Very Large Array (VLA) to search for H I 
21-cm line emission in the circumstellar environment of 6 Cephei. 

H I observations have been used previously to trace extended circumstellar emission sur- 



1 For the present work, we have recomputed the space motion vector for S Cephei using the updated solar 
constants from Schonrich et al. 2010. 
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rounding a number of mass-losing AGB stars (e.g., Gerard & Le Bertre 2006; Matthews & 
Reid 2007; Libert et al. 2008; Matthews et al. 2008, 2011). Frequently, the Hi emission 
associated with AGB stars is highly extended (up to ~1 pc) and shows signatures of in- 
teraction with the ISM in the form of trailing Hi wakes, velocity gradients caused by ram 
pressure effects, and/or density enhancements that demarcate the interstellar-circumstellar 
interaction zone. Furthermore, several of the evolved stars detected in H I are among those 
with FIR (and/or FUV) emitting bow shocks, including Mira, R Cas, and IRC+10216, un- 
derscoring the complementarity of these tracers for probing the chemistry and kinematics 
of stellar outflows and their interaction with their environments. This paper represents the 
first extension of this approach to the study of Cepheids. 



2. VLA Observations 

Hi 21-cm line observations of 5 Cephei were obtained using the VLA of the National 
Radio Astronomy Observatory (NRAOjj during four observing sessions in 2009 (see Ta- 
ble 2). All data were obtained using the compact (D) configuration (0.035-1.0 km baselines), 
providing sensitivity to emission on scales of up to ~15'. The primary beam of the VLA at 
the observing frequency of 1420.5 MHz is ~ 31'. 

The VLA correlator was used in 1A mode with a 1.56 MHz bandpass, yielding 512 
spectral channels with 3.05 kHz (~0.64 km s _1 ) spacing in a single (right circular) polariza- 
tion. The band was centered at a velocity of —7.5 km s _1 relative to the local standard of 
rest (LSR). In total, ~15.1 hours of integration were obtained on 5 Cephei. Approximately 
20% of the observed visibilities were flagged because of radio frequency interference (RFI) or 
hardware problems. Observations of S Cephei were interspersed with observations of a phase 
calibrator, 2203+626, approximately every 20 minutes. 3C48 (0137+331) was used as a flux 
calibrator, and an additional strong point source (2253+161) was observed as a bandpass 
calibrator (Table 2). To insure that the absolute flux scale and bandpass calibration were 
not corrupted by Galactic emission in the band, the flux and bandpass calibrators were each 
observed twice, with frequency shifts of +1.1 MHz and —1.7 MHz, respectively, relative to 
the band center used for the observations of 5 Cephei and 2203+626. 2203+626 was also 
observed once at each of these offset frequencies to permit more accurate bootstrapping of 
the absolute flux scale to the 5 Cephei data, although in some cases, one of the offset scans 
had to be discarded (see Table 2). We estimate that the resulting absolute flux scale has an 



2 The National Radio Astronomy Observatory is operated by Associated Universities. Inc., under cooper- 
ative agreement with the National Science Foundation. 
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Table 1. Coordinates and Stellar Properties of 5 Cephei 



Parameter 


Value 


Ref. 


a (J2000.0) 


22 29 10.2 


1 


5 (J2000.0) 


+58 24 54.7 


1 


I 


105?19 


1 


b 


+0?53 


1 


Distance (pc) 


273±0.011 


2 


Spectral Type 


F51b-Gllb 


3 


Pulsation period (days) 


5.366341 


3 


Mean T aB (K) 


5910 


4 


Mass (pulsation) 6 (M Q ) 


4.5±0.3 


5 


Mass (evolutionary) c (Mq) 


5.7±0.5 


5 


Mean My 


-3.47±0.10 


2 


Mean Luminosity (Z/q) 


~2000 




Mean radius (Rq) 


44.5 a 


6 


V d 

LSR 


-4.7 km s _1 


7 



a Based on the mean, limb darkened angular diame- 
ter 0ld=1-52O mas and the distance adopted here. 

b Based on V — K colors. 

c Assuming a canonical model with no overshoot; see 
§H3J. 

d Derived from the heliocentric radial velocity = 
-16.8 km s" 1 . 

Note. — Units of right ascension are hours, min- 
utes, and seconds, and units of declination are degrees, 
arcminutes, and arcseconds. All quantities have been 
scaled to the distance adopted in this paper. 

References. — (1) SIMBAD database; (2) Benedict 
et al. 2002; (3) Samus et al. 2011; (4) Andrievsky et 
al. 2005; (5) Caputo et al. 2005; (6) Armstrong et al. 
2001; (7) Wilson 1953 
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uncertainty of ~10%. 

The data were calibrated and imaged using the Astronomical Image Processing System 
(AIPS). At the time of our observations, the VLA contained 22 operational antennas with 
L-band receivers, 20 of which had been retrofitted as part of the Expanded Very Large Array 
(EVLA) upgrade. Data obtained during this EVLA transition period require special care 
during calibration^ After applying the latest available corrections to the antenna positions 
and performing an initial excision of corrupted data, we computed and applied a bandpass 
calibration to our spectral line data to remove closure errors on VLA-EVLA baselines. The 
bandpass was normalized using channels 163-448, thus excluding the portion of the band 
affected by aliasing. We next computed a frequency-averaged (channel 0) data set for use in 
calibrating the frequency-independent complex gains, again using channels 163-448. Follow- 
ing gain calibration, we applied time-dependent frequency shifts to the data to compensate 
for changes caused by the Earth's motion during the course of the observations. At this 
stage, we also applied hanning smoothing in velocity and discarded every other channel, 
yielding a 256 channel data set with a velocity resolution of ~1.3 km s _1 . 

Prior to imaging the line data, the u-v data were continuum-subtracted using a zeroth 
order fit to the real and imaginary components of the visibilities. Channels 20-95 of the 
hanning-smoothed data set (corresponding to an LSR velocity range of 35.6 to 132.2 km s _1 ) 
were determined to be line-free and were used for these fits. Although the spectral shape of 
the aliased portion of the continuum was better approximated by a higher order polynomial, 
the weakness of the continuum in the line data coupled with the lack of line-free channels on 
the high-frequency end of the band did not provide adequate constraints for a higher order 
fit. 

We imaged the H I line data using the standard AIPS CLEAN deconvolution algorithm 
and produced data cubes using various weighting schemes. Some characteristics of these 
cubes are summarized in Table 3. Additionally, we produced an image of the 21-cm con- 
tinuum emission in the 5 Cephei field using the line-free portion of the band. The peak 
continuum flux density within the primary beam was ~0.5 Jy (after correction for beam 
attenuation). We compared the measured flux densities of several sources in the field with 
those from the NRAO VLA Sky Survey (Condon et al. 1998) and found agreement to 
within <7%. We did not detect any continuum emission at or near the position of 5 Cephei 
and place a 3a upper limit on the 21-cm continuum flux density at the stellar position of 
Fcont <1.0 mJy This is consistent with the results of Welch & Duric (1988), who found no 
evidence of significant mass-loss from S Cephei in the form of an ionized wind based on radio 



http://www.vla.nrao.edu/astro/guides/evlareturn/ 
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continuum observations at 5 GHz. 



3. Results 

Figure Q] shows a series of H I channel maps extracted from our Tapered VLA data cube 
(see Table 3). Only a few representative channels are shown blueward of the stellar systemic 
velocity (Vlsr = —4.7 km s^ 1 ), as this portion of our band was strongly contaminated by 
large-scale Galactic emission along the line-of-sight. This emission is poorly spatially sampled 
by the VLA, resulting in patterns of positive and negative mottling across the channel images. 
To better illustrate the character of the Galactic emission, Figure [2] presents a single-dish 
Hi survey spectrum from Kalberla et al. (2005), which shows Hi emission with a brightness 
temperature as high as ~85 K toward the direction of 5 Cephei. Both the strength and 
spatial complexity of this emission hampered our search for circumstellar signals over the 
velocity range —125 <Vlsr <12.5 km s _1 . 

Fortunately, the spectrum in Figure [2] reveals that the Galactic H I emission drops pre- 
cipitously for Vlsr >13 km s _1 , providing a clean band over which to search for circumstellar 
emission. Because the outflow velocities from Cepheids are expected to be at least a few 
tens of km s _1 (e.g., Deutsch 1960; Reimers 1977; Holzer & MacGregor 1985), this velocity 
window is of prime interest for searching for circumstellar emission associated with 5 Cephei. 

Consistent with the total power spectrum in Figure [21 we see in Figure [1] that the 
pattern in the VLA images caused by large-scale Galactic emission diminishes near Vlsr — 
13.8 km s -1 . Furthermore, that same channel exhibits an extended Hi emission region 
centered near the position of 6 Cephei. Looking toward higher velocities, we find a series 
of 15 contiguous channels where spatially extended H I emission is detected coincident with 
the position of 5 Cephei. Several marginally extended features of ~ 3a significance are also 
seen clustering near the stellar position in each of the next five contiguous channels, out to 
VisR = 38.2 km s" 1 . 

In Figure [3] we present velocity-integrated Hi total intensity contours, derived from data 
spanning the velocity range 13.8 < Vlsr < 38.2 km s _1 . These contours are overlaid on the 
Spitzer 24yum and 70/im maps, respectively, from Marengo et al. (2010b). Some caution is 
required in interpreting the detailed morphology of the velocity-integrated Hi distribution 
in Figure [31 as additional material associated with this nebula is likely to be present at lower 
velocities, where it cannot be disentangled from the prominent Galactic emission. However, 
despite this caveat, several characteristics of the integrated Hi emission seen in Figure [3] are 
strongly indicative of an association with the circumstellar environment of S Cephei. 
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Table 2. VLA Calibration Sources 



Source 


a(J2000.0) 


<5(J2000.0) 


Flux Density (Jy) 


Date 


3C48 a 


01 37 41.29 


+33 09 35.13 


15.88* 


All 


2253+161 b 


22 53 57.74 


+16 08 53.56 


13.74±0.15t 


2009-Oct-28 








13.46+0.28** 


2009-Nov-15&16 








13.27+0.23** 


2009-Nov-22&23 








13.36±0.22t 


2009-Nov-27&28 


2203+626 c 


22 03 20.96 


+62 40 34.27 


2.75+0.02" 


2009-Oct-28 








2.72+0.01" 


2009-Nov-15&16 








2.73+0.01" 


2009-Nov-22&23 








2.75±0.0l t1 " 


2009-Nov-27&28 



Note. — Units of right ascension are hours, minutes, and seconds, and units of 
declination are degrees, arcminutes, and arcseconds. 

'Adopted flux density at 1420.5 MHz, computed according to the VLA Calibration 
Manual (Perley & Taylor 2003). 

t Computed flux density at 1419.2 MHz; data at 1422.1 MHz were corrupted. 

"Mean computed flux density from observations at 1419.2 and 1422.1 MHz; see §[2] 

ttComputed flux density at 1422.1 MHz; data at 1419.2 MHz were corrupted. 

a Primary flux calibrator 

b Bandpass calibrator 

c Phase calibrator 



Table 3. Deconvolved Image Characteristics 



Image 


n 


Taper 


#FWHM 


PA 


rms 


Descriptor 




(kA,kA) 


(arcsec) 


(degrees) 


(mJy beam -1 ) 


(1) 


(2) 


(3) 


* (4) ' 


(5) 


(6) 


Robust +1 


+1 




49" x 45" 


20°3 


1.6 


Natural 


+5 




54" x 49" 


21°7 


1.5 


Tapered 


+5 


2,2 


97" x 89" 


2°5 


1.8 


Continuum 


-1 




39" x 36" 


0°3 


0.39 



Note. — Explanation of columns: (1) image or data cube designation 
used in the text; (2) AIPS robustness parameter used in image deconvo- 
lution; (3) Gaussian taper applied in the u and v directions, expressed as 
distance to 30% point of Gaussian in units of kilolambda; (4) dimensions of 
synthesized beam; (5) position angle of synthesized beam (measured east 
from north); (6) mean rms noise per channel (1 a) in the unaliased portion 
of the band (line data) or in frequency-averaged data (continuum) . 
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One of the most intriguing aspects of the H I emission distribution in Figure [3] is its 
"head-tail" morphology. This is reminiscent of the structures of the extended circumstellar 
envelopes of several of the AGB stars that have been imaged in Hi, FIR and/or FUV 
light (Martin et al. 2007; Matthews et al. 2008, 2011; Ueta 2008; Sahai & Chronopoulos 
2010; Jorissen et al. 2011) and results when material shed by a stellar wind is swept by ram 
pressure into a trailing wake as the star transverses the ISM (e.g., Wareing et al. 2007b). The 
direction of space motion of 5 Cephei, indicated by an arrow on Figure[3], strongly supports an 
analogous interpretation for the H I nebula surrounding this star. We see that the maximum 
angular extent of the H I nebula (~ 13', corresponding to ~1.0 pc at the distance of 5 Cephei) 
is aligned with the space motion direction. Furthermore, the distribution of the H I nearest 
to the position of 5 Cephei (r <2') is well correlated with the distributions of both the 24/zm 
and 70/zm emissions detected by Spitzer. The Hi emission thus appears to be associated 
with the same dusty wind traced by the IR data. 

The peak of the Hi column density in Figure |3] (iVm ~ 2.0 x 10 19 cm~ 3 ) does not 
correspond to the position of 6 Cephei or its companion HD 213307 (indicated on Figure [3] 
by a star and cross symbol, respectively), but rather lies ~ 112" to the northeast of S Cephei, 
near the apex of the shock structure delineated by the IR emission. Such an offset of the 
peak emission from the star is not surprising, as this kind of leading density enhancement 
is a hallmark of the "snowplow" effect that may occur when a wind-emitting star moves 
supersonically through the ambient medium (Isaacman 1979; Raga et al. 2008). 

To the west of 5 Cephei, i.e., downstream from the star's trajectory, our VLA image 
reveals a wake of material trailing the motion of the star. In contrast to the gaseous wakes 
associated with stars such as Mira (Matthews et al. 2008) and X Her (Matthews et al. 
2011), the material downstream from <5 Cephei is not well collimated and more closely re- 
sembles the wide-angle "vortical tail" associated with the carbon star IRC+10216 (Sahai & 
Chronopolous 2010). Numerical simulations of the interactions between stellar winds and 
streaming environments frequently predict broad, turbulent wakes, qualitatively similar to 
what we see in Figure [3] (e.g., Wareing et al. 2007a), although the degree of collimation and 
other properties of circumstellar wakes are highly sensitive to the stellar velocity, properties 
of the stellar wind, and the density of the ambient medium (Comeron & Kaper 1998), such 
that fully unraveling the interplay between these factors for a particular source generally 
requires a specifically tailored model. 

To gauge the total mass of H I giving rise to the emission seen in Figures [1] and [3] over 
the velocity range Vlsr=15.0 to 38.2 km s -1 , we have summed the emission surrounding 
the stellar position (after correction for primary beam attenuation) in irregularly shaped 
blotches defined by 3cx brightness contours. This approach yielded an integrated Hi flux 
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density of / S v dv = 4.1 ± 0.2 Jy km s _1 . Using the relation M ffl = 2.36 x 10" 7 d 2 / S v dv M Q 
where d is the stellar distance in pc (e.g., Roberts 1975), this translates to an Hi mass 
of Mhi ~0.07M q . While this approach for measuring the total Hi context minimizes the 
contribution of background noise to the measurement, it may result in the exclusion of 
weak emission or emission from regions not contiguous with the primary emission region in 
each channel. For comparison, we therefore summed the emission in each channel over a 
fixed rectangular aperture (16'6xl3'3, centered 100" west of 5 Cephei's position). We found 
S v = 4.4 ± 0.4 Jy km s _1 , consistent with our first estimate. 

We interpret our measured integrated H I flux density as a lower limit to the amount of 
circumstellar material present around 5 Cephei for two reasons. First, additional emission 
may lie undetected at velocities where our spectral images were contaminated by Galactic 
emission. For example, if we assume that the H I line profile is symmetric about the systemic 
velocity of 5 Cephei, this would imply that roughly two-thirds of the channels containing 
circumstellar signal were excluded from our flux density measurement. Secondly, at our 
observing frequency, the VLA D configuration is only sensitive to emission on angular scales 
<15'. Since this is comparable to the angular extent of the observed Hi nebula surrounding 
S Cephei, it is plausible that we have missed additional extended emission. This problem 
could be exacerbated if the most spatially extended material is decelerated as a result of its 
interaction with the ISM (e.g., Matthews et al. 2008); this could in turn cause its velocity 
to be blueshifted into the portion of our band contaminated by Galactic emission. Mapping 
the region in H I with a single-dish telescope would be of considerable interest as a means of 
searching for an additional large-scale emission component. 



4. Discussion 

4.1. How Robust is the Association between the Hi Nebula and 5 Cephei? 

The Galactic ISM is well known to exhibit structures on wide variety of scales, down to 
arcminutes and below (e.g., Verschuur 1974; Greisen & Liszt 1986; Heiles & Troland 2003; 
Begum et al. 2010). We therefore have considered the possibility that the Hi nebula that 
we have detected toward 5 Cephei could simply be a chance superposition of an unrelated 
Hi cloud along the line-of-sight. However, several lines of evidence indicate that this is 
extremely unlikely. 

As described above, the close correspondence between the Hi emission contours sur- 
rounding 5 Cephei and the bow shock structure traced by FIR emission implies that the two 
are linked. Since the formation of a bow shock requires a stellar wind, a linkage between the 
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H I and FIR emission argues strongly against the former arising from a random interstellar 
cloud. 

Based on a survey with the Arecibo L-band Feed Array, Begum et al. (2010) identified 
a population of compact Hi clouds associated with the Galactic plane (but distinct from 
the main Galactic disk emission) whose column densities and angular sizes overlap with 
those characterizing the Hi nebula that we have detected toward S Cephei. However, the 
mean characteristic linewidth of the clouds catalogued by Begum et al. is AV—A.2 km s _1 
(FWHM), with values generally lying in the range of ~l-8 km s -1 . Only a handful of the 
96 clouds in their sample exhibit velocity widths comparable to or larger than that of the 
integrated linewidth of the 5 Cephei nebula (Ay 12.6 km s _1 based on a Gaussian fit). 
Furthermore, while Begum et al. have underscored the possibility that a subset of their 
compact clouds may in fact be associated with mass-loss from variable stars, the closest 
association that they identified between a catalogued variable and a cloud from their sam- 
ple has a projected separation of 10'. Since the Hi properties of the 8 Cephei nebula are 
somewhat atypical of the general population of compact Galactic clouds, we conclude that 
the probability of detecting such a cloud with a head-tail morphology at both the position 
and velocity expected for S Cephei, is extraordinarily small. 



4.2. The Stellar Mass-Loss Rate 

4-2.1. New Estimates of M Based on the VLA Data 

One of the most powerful applications of our VLA measurements is that they allow 
us to obtain new, independent constraints on the mass-loss rate from S Cephei. Two key 
advantages of Hi observations for characterizing stellar mass-loss are first, that the Hi 
mass can be derived directly from the observed emission (assuming it is optically thin) and 
secondly, that spectroscopic imaging of Hi is sensitive to material much farther from the 
star than most other observational techniques for probing stellar mass loss, thereby tracing 
a larger fraction of the total mass-loss history. 

We consider two independent approaches for estimating the mass loss rate of 5 Cephei 
from our VLA data. Our first method is to adopt a simplified model in which the star is 
assumed to be undergoing mass-loss through an isothermal, spherically symmetric wind that 
is in free expansion out to some radius r = Rq. While pulsation-driven winds are inherently 
non-steady, we assume that the outflow can be approximated as quasi-steady when averaged 
over the timescales sampled by our observations. 

In the case where the H I emission is optically thin, a theoretical line profile as a function 
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of velocity, V, for this idealized wind can be written as: 



S(V) 



2/Vvrln 2 M 
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(see Knapp & Bowers 1983; Olofsson et al. 1993). Here / is the conversion factor between 
units of brightness temperature in kelvins and mJy per beam for the radio telescope used, d 
is the distance to the star, B is the half-power synthesized beamwidth of the telescope, V Q is 
the outflow velocity of the wind, /x is the mean atomic weight of the gas (taken to be 1.3), 
and tuh is the mass of the hydrogen atom. All quantities are expressed in cgs units. 

For comparison with the theoretical profile, we extracted from the VLA data a spectrum, 
averaged over a single synthesized beam centered on 5 Cephei. We used the "Robust+1" data 
cube for this analysis (Table 3), which has the highest angular resolution and should minimize 
contamination from structures in the H I envelope that are expected to lie outside the freely 
expanding wind (e.g., material in the tail and/or bow shock). We then compared the observed 
spectrum with a series of models where the mass-loss rate and wind outflow speed were taken 
as free parameters. Before comparison with the data, each model computed from Eq. 1 was 
convolved with a broadening kernel to account for the finite spectral resolution and the gas 
turbulence. The model parameters were constrained by scanning a grid of possible values and 
minimizing x 2 - The outflow velocity, mass-loss rate, and smoothing kernel were incremented 
by values of 0.1 km s _1 , 1 x 10 -8 M yr^ 1 , and 0.1 km s _1 , respectively. The best-fitting 
model profile is overp lotted on the observed spectrum in Figure HI The model parameters are 
V = 35.6 ±1.2 km s~\ M = (1.7 ±0.1) x 10~ 6 M yr" 1 , and a turb = 1.8 ±0.2 km s _1 . The 
quoted uncertainties do not reflect the uncertainties in parameters held fixed in the model, 
namely the stellar systemic velocity, the envelope radius, and the distance to the star. Our 
model underestimates the peak H I flux density of the observed spectrum, which may reflect 
uncertainties in the fixed parameters of the model (e.g., the envelope radius) or be due to 
contamination of the observed spectrum by swept-up interstellar gas and/or gas from outside 
the free expansion zone of the wind. It may also reflect a breakdown of our highly idealized 
model (e.g., owing to highly episodic mass-loss or a latitude dependent wind). 



A second, independent means of constraining the mass-loss rate from the H I data comes 
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from the total mass of circumstellar H I emission observed, coupled with available constraints 
on the timescale for the mass-loss. Over the velocity range 13.8 < Vlsr < 38.2 km s _1 , we 
directly measured a total Hi mass of 0.07 M from material that appears to be associated 
with the circumstellar environment of S Cephei (see §[3]). Applying a multiplicative correction 
of 1.34 to account for the mass of helium yields M tot >0.09 M Q . As noted above, we regard 
this as a lower limit to the total envelope mass, since additional circumstellar emission may 
lie in portions of our band that are contaminated by Galactic signal and/or at large angular 
scales to which the VLA is insensitive. 

The tail of H I emission downstream from 5 Cephei has an extent of ~ 500", or 2 x 10 18 cm 
in the plane of the sky. Assuming a tangential velocity of 7 km s -1 (derived from the 
observed stellar radial velocity and proper motions as described in Marengo et al. 2010b 
and the new solar constants of Schonrich et al. 2010), the crossing time for the nebula is 
~90,600 yr. For the total circumstellar mass above, this would imply a mean mass-loss rate 
of M ?a 1 x 10 -6 M Q yr -1 . While this is comparable to the rate derived from our model 
fitting above, the actual mass-loss rate could be smaller if some fraction of the observed H I 
gas arises from swept-up interstellar material rather than directly from stellar mass-loss (see 
§ l4.2.2l below). Moreover, the dynamical crossing time of a stellar wake provides only a lower 
limit to its age (see Matthews et al. 2008); if 5 Cephei has been losing mass over a more 
prolonged interval, M would also be lower. We note, however, that 6 Cephei is believed to 
be on its second crossing of the Cepheid instability strip (Berdnikov et al. 2000; Turner et al. 
2006), and the total duration of this phase for a star of its mass and metallicity is expected 
to be ~(1.3-5.9) xlO 5 yr (see § 14.51) . These numbers are only a few times larger than our 
estimated dynamical age, suggesting that the mean mass loss rate cannot be significantly 
lower than we have inferred under the assumption that all of the observed mass-loss occurred 
during the Cepheid phase. [The duration of the first Cepheid crossing (~ 15, 000 — 30, 000 yr 
for a star like 5 Cephei; see § 14.51) is too short to impact this conclusion]. On the other hand, 
if we have underestimated the amount of emission in the Hi nebula because of Galactic 
confusion, this would imply a higher mass-loss rate, although this correction is expected to 
be less than a factor of ~2 or 3, assuming the line profile is symmetric. 

4-2.2. Is a Significant Fraction of the Circumstellar Material Swept from the ISM? 

As mass-losing stars move through the ISM, their jetsam can be augmented by material 
swept from the local ISM. The fraction of the circumstellar envelope that originates in this 
manner will depend on a combination of the stellar mass, the stellar wind parameters, and 
the ambient ISM density (Villaver et al. 2002). Accounting for this swept-up interstellar 
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mass is important in evaluating the true mass-loss rate from the star. 

In the absence of a detailed numerical model, we can obtain a crude estimate of 
the amount of swept-up material in the circumstellar environment of 5 Cephei as M s = 
J ' m ' 7rr 2 V s p dt, where r is the effective radius of the star, V s is its space velocity, p is the 
mass density of the ambient medium, and t m i is the length of time that the star has been 
losing mass. 

Adopting an ISM particle density n ~10 cm -3 (see § 14.2. 3p translates to p = 1.7 x 
1CT 23 g cm -3 . Assuming r=0.107 pc (Marengo et al. 2010b), V,=10.3 km s" 1 (§ HJ, and 
t mi =90,600 yr (§ IP]) , this yields M s w 0.01(ra ffl /10 cm~ 3 )(t m/ /90, 600 yr) M , or roughly 
14% of the directly observed H I mass. Even after allowing for factors of a few uncertainty 
in both n and t m[ , this supports our interpretation of the emission detected by the VLA as 
arising predominantly from a wind from 5 Cephei. This is also consistent the conclusions 
of Marengo et al. (2010b), who argued that the low inferred PAH content of the material 
detected by Spitzer implies that it is dominated by debris originating in the stellar wind. In 
the future, more sophisticated numerical simulations (cf. Villaver et al. 2002) could help to 
obtain more accurate constraints on M s . 



4-2.3. Comparison with Previously Derived Mass-Loss Rates for 5 Cephei 

Based on their Spitzer measurements, Marengo et al. (2010b) used two indepen- 
dent approaches to estimate a mass-loss rate for 5 Cephei in the range M ~ 5 x 10~ 9 
to 6 x 10~ 8 M Q yr -1 . If we scale these values to account for revised stellar space ve- 
locity of V s =10.3 km s _1 adopted in the present work (see § [T]) and the outflow velocity 
of V Q ~35.6 km s _1 inferred from the Hi data, this translates to M w 7.2 x 10 -9 and 
2.1 x 10~ 8 M© yr~ x , respectively. These values are consistent with theoretically predicted 
mass- loss rates for Cepheids (Neilson & Lester 2008), but are significantly lower than the 
values that we infer from the H I data. 

The smaller of the two mass-loss rate values from Marengo et al. (2010b) was derived 
based on the observed stand-off distance of the bow shock structure together with ram 
pressure balance arguments. The implied mass-loss rate could be higher if the local ISM 
particle density in the vicinity of 5 Cephei is larger than the originally adopted value of 
n =0.55 cm -3 . Indeed, there is evidence that a higher particle density in the region is 
plausible. Fernie (1990) derived a reddening value of E(B — V) — 0.092 toward 5 Cephei; 
assuming AyjE[B — V)=3.1 (Savage & Mathis 1979), this implies a \^-band extinction 
coefficient Ay w 0.29. Such a value is typical of interstellar clouds with n ^10-50 cm -3 
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(Turner 1997). We note that owing to its larger angular extent, the circumstellar dust 
detected by Spitzer should contribute negligibly to Ay. Assuming a typical grain radius 
a=0.001yum, a total dust mass M dust =6 x 1CT 7 M , and a radial extent r=0.103 pc (Marengo 
et al. 2010b), then Ay^ oca x ~ 0.22Md ust a _1 r -2 0.01. Finally, from the database of Kalberla 
et al. (2005), the Hi column density along the line-of-sight to 8 Cephei is Nm = 7.7 x 
10 21 cm -2 , implying a mean particle density out to the distance of the star of n ^9.2 cm -3 . 
This is again consistent with a local particle density several times higher than the canonical 
Galactic plane value for warm, intercloud Hi (cf. Dickey & Lockman 1990). Assuming a 
local ISM density with n ~ 10 cm~ 3 would imply M ~ 1.4 x 10~ 7 M Q yr _1 based on ram 
pressure balance arguments, in somewhat better agreement with our H i-derived value given 
the uncertainties in both estimates. 

The second mass-loss rate estimate from Marengo et al. (2010b) was derived under 
the assumption that the gas-to-dust (g/d) ratio of the circumstellar material has a value 
of 100, typical of what is found in the ISM and in the circumstellar environments of AGB 
stars. However, it is presently unknown whether this value should hold in the circumstellar 
environment of a Cepheid. For example, the ratio g/d is found to be 5-10 times higher 
than the canonical value in the circumstellar envelopes of some supergiants (e.g., Skinner & 
Whitmore 1988), presumably because grain condensation has been incomplete in the winds 
of such stars (Mauron & Josselin 2011). Even more extreme values may be plausible in the 
winds of Cepheids owing to their higher temperatures and higher ultraviolet fluxes. In the 
case of 5 Cephei, g/d « 2300 would be needed to reconcile the mass-loss rate of Marengo et 
al. (2010b) with our Hi estimates. Interestingly, such a value would be self-consistent in the 
sense that the empirically derived mass volume density of the dusty wind, pof, would then 
match the density for a wind with M 1 x 10~ 6 M© yr _1 and V Q =35.6 km s _1 as predicted 
by the mass continuity equation for stellar winds: M = 47ipopr 2 V (Lamers & Cassinelli 
1999). 

In summary, our best estimate of the mass-loss rate of 5 Cephei based on the combined 
results from our VLA Hi data and Spitzer infrared data is M m (1 ± 0.8) x 10 -6 M Q yr _1 . 
The quoted error bar reflects the dominant systematic uncertainties and is based on the 
dispersion in the values derived from the independent methods presented in the preceding 
sections. 
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4.3. Implications of the Empirically Derived Mass-Loss Rate for Theoretical 

Models of Cepheids 

4-3.1. Consistency with Theoretically Predicted Mass Loss Rates 

The majority of Cepheids for which mass-loss rates (or upper limits) have previously 
been determined (§[[]) are inferred to have M values at least two or three orders of magnitude 
smaller than we have derived for 6 Cephei in the present work. We emphasize however, 
that none of these previous determinations exclude the possibility of significant mass-loss 
through a predominantly neutral atomic wind, of the type we have detected through Hi 
observations. In the case of 5 Cephei, the wind also revealed its presence through FIR 
emission, although this need not be the case for all mass-losing Cepheids, particularly those 
with higher effective temperatures (where dust grains might be destroyed more readily), those 
that lack a companion (which may play a role in exciting FIR emission), or those residing in 
lower density environments (where dust swept from the ISM might contribute negligibly to 
the circumstellar FIR emission) . It is clear that an H I survey of a larger sample of Cepheids 
will be needed to obtain a better understanding of the importance of these factors and to 
assess the ubiquity of Cepheid mass loss that is traceable through Hi 21-cm line emission. 

Based on theoretical studies to date, the most promising physical mechanism by which 
Cepheids could lose significant mass appears to be the enhancement of radiatively driven 
mass loss through shocks and pulsation (Willson & Bowen 1984; Neilson & Lester 2008). 
From analytic models of this type, Neilson & Lester (2008) and Neilson et al. (2011) predicted 
typical Cepheid mass-loss rates in the range 10~ 10 to 10 -7 M Q yr _1 , although these rates 
are highly sensitive to the balance of forces and the adopted physical parameters for a 
given star. In particular, for shorter periods (smaller masses), the mass-loss rate has a non- 
linear dependence on the stellar mass. Thus for a Cepheid with the mass and period of 
S Cephei, it appears that the models of Neilson & Lester and Neilson et al. can in principle 
accommodate mass-loss rates as high as 10 -6 M & yr _1 , consistent with our VLA findings 
(H. Neilson, private communication). However, this conclusion may change if mass loss is 
highly episodic (see below). 

4-3.2. The Effects of Mass Loss on the Pulsation Period 

Secular period changes are well known to occur in Cepheids, but are generally assumed 
to be dominated by evolutionary effects (which result in changes in the effective temperature 
and luminosity) rather than mass loss (Turner et al. 2006). It is therefore of considerable 
interest to compare the change in pulsation period for 5 Cephei, P, that is predicted to 
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occur as a result of its derived mass loss rate with the value of P previously derived from 
observational data. 

Equation 62 of Neilson & Lester (2008) provides an analytic formula for estimating the 
effect of mass loss on the rate of change in the stellar pulsation period. The most extreme 
period changes due to mass loss are expected to occur when changes in both T e ff and L (and 
hence R) are minimal. Moreover, the time scale for the decrease in total mass is predicted 
to be faster than the canonical evolutionary time scale. Under these conditions, the period 
changes dominated by mass-loss should always be positive. If we take the pulsation period 
and pulsation mass for 8 Cephei from Table 1, together with M = 10 -6 M Q yr _1 , using 
Equation 62 of Neilson & Lester we estimate the period change as a result of mass loss 
from 5 Cephei to be P = 0.06 s yr^ 1 . This is in apparent contradiction to the observed 
value of P — —0.1 s yr _1 derived for 5 Cephei by Berdnikov et al. (2000) based on data 
spanning 130 years. This suggests that mass loss is only one contributor to the secular period 
change and that evolutionary effects on the star following mass loss cannot be ignored. An 
additional consideration is that mass loss inside the instability strip may be episodic. This 
could explain, for example, why most previous searches for Cepheid mass loss (which have 
only been sensitive to material in close proximity of the star) have typically implied much 
lower rates of mass-loss than we infer from the VLA and Spitzer observations of S Cephei. 
However, in order to account for the same total amount of matter shed, highly episodic mass 
loss would require a boost in efficiency of the Cepheid mass loss mechanism. 



4-3.3. Could the Mass Loss Have Occurred During an Earlier Evolutionary Phase? 

One of the few Cepheids that has been previously inferred to have a mass loss rate 
comparable to the rate we have derived for 5 Cephei in our present study is RS Pup, where 
M ps 10~ 6 M Q yr _1 (Deasy 1988; Barmby et al. 2011). RS Pup is unique among Galactic 
Cepheids in that it is surrounded by a prominent reflection nebula. The linear extent of 
this nebula (~1 pc) is comparable to the Hi nebula that we measure toward 5 Cephei. 
Although the origin of RS Pup's nebulosity is consistent with copious mass loss, some authors 
have argued that the bulk of the nebular material may have originated during an earlier 
evolutionary phase, either as a Be star (Kervella et al. 2009) or as a red giant (Havlen 1972). 
The question of whether such arguments could also be applied to 5 Cephei is important, 
since it impacts whether subsequent mass loss during the Cepheid phase will be sufficient to 
reconcile the Cepheid mass discrepancy for this star. 

The geometry of the nebula that we observe surrounding 5 Cephei places strong con- 
straints on this issue. The presence of an infrared-emitting bow shock and a head-tail struc- 
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ture to the H I nebula strongly suggest that 5 Cephei has ongoing mass-loss. While this alone 
does not exclude the possibility that some portion of the material in its circumstellar nebula 
was shed during an earlier evolutionary phase (e.g., at the tip of the red giant branch), the 
fact that the dynamical age of the H I nebula is smaller than the predicted duration of the 
second crossing of the instability strip for 5 Cephei (see § 14. 5[) suggests a strong likelihood 
that the bulk of the mass loss occurred during the Cepheid evolutionary phase. Further, the 
dynamical age of the nebula implies that even if some of the mass were lost during the red 
giant phase, this is likely to have immediately preceded the current Cepheid crossing. Hence 
the resulting mass loss would still necessitate a change in the evolutionary mass of the star 
between its first and its second (current) Cepheid crossing. 

4.4. Comments on the Inferred Outflow Velocity for the S Cephei Wind and 

Its Implication for the Mass Loss Mechanism 

In § 14.2. 1| we estimated an asymptotic (outflow) velocity for the 6 Cephei wind of 
V ~35.6 km s -1 . To our knowledge, this is the first directly measured outflow velocity 
for a Cepheid. An important characteristic of this value is that it is significantly smaller 
than the stellar escape velocity. Adopting the mean stellar radius and pulsation mass for 
5 Cephei quoted in Table 1, the escape velocity from the photospheric surface of 5 Cephei is 
V^, sc ~ 200 km s -1 , or several times the asymptotic velocity. While the radius and effective 
surface gravity of 5 Cephei change over the course of its pulsation cycle, these effects are 
not sufficient to change this conclusion. Andrievsky et al. (2005) reported changes in the 
effective surface gravity of 5 Cephei by a factor of 5.6, which implies that even when the 
surface gravity is at its minimum, the escape velocity would still be V^ sc ~ 120 km s _1 . 
S Cephei thus appears to share one of the defining characteristics of the cool winds of AGB 
stars and late-type supergiants, where measured outflow velocities are in general several 
times smaller than the photospheric escape speed (i.e., V « V esc ; see Holzer & MacGregor 
1985; Judge & Stencel 1991). Interestingly, 5 Cephei also appears to adhere to the empirical 
relation V Q ~ 1.6 x 10~ 3 V^ C found by Reimers (1977) based on observations of K giants 
and G and K supergiants (see also Judge 1992), although the significance of this is unclear, 
since the physical underpinnings of this relation are poorly understood (Reimers 1977) and 
5 Cephei is not expected to be undergoing mass loss via the same mechanism as the stars in 
Reimers' sample. 

As discussed by Holzer & MacGregor (1985), the condition that V Q « V esc has im- 
portant implications for understanding and constraining the mechanism(s) driving a stellar 
wind. For example, this requires the existence of some type of regulatory processes such that 
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most of the driving energy of the wind goes into lifting material out of the stellar gravita- 
tional field rather than accelerating the flow. It follows that the inferred momentum flux of 
the wind, MV a , represents only a small fraction of the overall energy required to drive the 
wind. An additional implication is that most of the energy added to the wind must be in 
the form of momentum rather than heat. 



4.5. A Possible Solution to the Cepheid Mass Discrepancy for 5 Cephei? 

A fundamental question raised by our new empirically derived mass-loss rate for S Cephei 
is whether this degree of mass loss can reconcile the mass discrepancy for the star. As 
described above, 8 Cephei is believed to be on its second crossing of the instability strip. 
The total duration of this evolutionary stage (as well as for the first and third crossings) has 
previously been calculated for Cepheids with different masses and metallicities by Bono et 
al. (2000a). However, these models did not include the effects of mass loss and therefore 
their direct application of the case of 5 Cephei will not yield self-consistent results. To 
overcome this problem, we computed new intermediate-mass evolutionary models specifically 
for S Cephei. The theoretical framework we employed was developed by Pietrinferni et al. 
(2004). The evolutionary tracks we adopted are available in the BaSTI data bas^l and 
were constructed by assuming a scaled solar chemical mixture and a fixed solar chemical 
composition with helium and metal abundances of F=0.273 and Z=0.02, respectively. 

Several lines of empirical evidence support the occurrence either of mild convective 
core overshooting during the central hydrogen burning phases of intermediate mass stars 
(Barmina et al. 2002; Cordier et al. 2002; Cassis & Salaris 2011; Evans et al. 2011) or of 
enhanced post-main sequence mass loss (Bono et al. 2002; Natale et al. 2008; Cantiello et 
al. 2011). However, to constrain the effects of different mass loss rates, our current canonical 
evolutionary models were constructed by neglecting extra-mixing sources both at the edge 
of the convective core and at the bottom of the convective envelope (Alongi et al. 1991). We 
plan to provide a more detailed discussion of the impact of the mass loss on the evolutionary 
properties of helium burning intermediate-mass stars in a forthcoming paper (Bono et al., 
in preparation). 

Figure [5] shows H-R diagrams for evolutionary tracks at a fixed initial stellar mass 
(M* = 5M ) but different assumptions concerning the mass-loss rate. Each panel shows the 
evolutionary phases just before and soon after helium burning, i.e. the so-called "blue loop" . 
The solid lines denote the predicted first overtone blue (hot) edge and the fundamental red 



4 A few selected evolutionary tracks were specifically computed by A. Pietrinferni for this project. 
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(cool) edge of the Cepheid instability strip according to Bono et al. (2000b). The triangle 
and the diamond mark the onset (tip of the red giant branch) and the end of central helium 
burning phases, respectively. The evolutionary time spent inside the instability strip is also 
indicated. The evolutionary track plotted in the top panel was constructed by neglecting 
mass loss. The track plotted in the middle panel was constructed by assuming a mass loss 
starting from the main sequence phase according to the semi-empirical formula from Reimers 
(1975): 

M = 4 x 10 -i^W!« (2) 
' M/M v ; 

The free parameter rj was fixed at 0.4 (Pietrinferni et al. 2004; Bono et al. 2006). Following 
this prescription, the mass loss rate during the evolution of the star will vary according 
to the changes in surface luminosity and radius. However, a comparison of the theoretical 
predictions plotted in the two top panels of Figure |5] shows that the impact of the mass loss 
based on the prescription of Reimers is negligible, and indeed, both the extent in effective 
temperature of the blue loop and the evolutionary time spent inside the instability strip are 
minimally affected. Plausible increases in the free parameter rj do not change this outcome, 
and the Reimers model cannot produce sufficient mass loss to account for the observed 
quantity of circumstellar material around 5 Cephei, nor can it reconcile the mass discrepancy 
for this star. 

Finally, in the bottom panel of Figure we show an evolutionary track with the same 
initial mass of the top panels, but assuming a steady mass loss rate of M = 5 x 10~ 7 M© yr _1 
during the evolutionary phases inside the Cepheid instability strip. As expected, we found 
that implementing such mass-loss causes a significant decrease in the extent in effective 
temperature of the blue loop owing to the decrease in the envelope mass (Bono et al. 2000a, 
and references therein), and in turn a decrease in the evolutionary time spent inside the 
instability strip. We also explored models with higher mass-loss rates (not shown), but found 
that for M >10 -6 M & yr" 1 , the temperature extent of the blue loop becomes significantly 
smaller, leading to implausibly short Cepheid lifetimes. 

To further constrain the impact of the initial stellar mass, we computed an additional set 
of evolutionary models with an initial stellar mass M* = 6M Q (Figure |6]). The evolutionary 
tracks show a behavior similar to the M* = 5M models when moving from the case of no 
mass-loss (top) to the Reimers mass-loss case (middle) to the case of mass loss inside the 
instability strip at a rate of M = 5 x 10 -7 M Q yr _1 (bottom). As with the M* = 5M 
models, higher mass-loss rates (M >10~ 6 M & yr _1 ) destroyed the blue loop, resulting in a 
negligible lifetime on the instability strip. 
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The above findings indicate that either (1) the true mass-loss rate from 5 Cephei is at 
the low end of our empirically derived range, or (2) that a more complex treatment of mass 
loss is required (e.g., episodic mass loss or implementation of multi-dimensional Eulerian 
hydrodynamic and stellar evolution codes; Mocak et al. 2009). The inclusion of extra 
mixing and stellar rotation may also be important. We consider possibility (1) below, while 
possibility (2) will be explored in future works. 

If we assume that 5 Cephei is currently near the midst of its second crossing of the 
instability strip, the models presented in Figures \5\ and E] predict that a Cepheid with the 
chemistry and evolutionary mass of 5 Cephei and undergoing mass loss at a rate of M — 
5 x 10~ ' M Q yr _1 will have so far spent roughly (0.9-3.0) x 10 5 yr on the instability strip (most 
of it on the second crossing). These numbers are consistent with the dynamical age of the 
H I wake estimated in § 14.2.11 Steady mass-loss at the prescribed rate over such an interval 
is predicted to result in a mass of circumstellar material of ~0.04 to O.15M . The upper end 
of this range is consistent with our empirically derived estimate for the mass of circumstellar 
material after correction for He (§ 14.2.11) . implying that such a model is self-consistent. 
Moreover, this amount of material is remarkably close to the change in mass required to 
reconcile the pulsation and evolutionary masses of 5 Cephei (M e — M p 0.2 — 1.8 M Q 
assuming that M e = 5.5 ± 0.5 M . We conclude that to within current observational and 
theoretical uncertainties, mass loss offers a tenable solution to resolving a significant fraction 
of the the mass discrepancy for 5 Cephei. 

From observations of only a single star, we cannot draw any general conclusion about the 
ubiquity or evolutionary importance of mass loss from Cepheids. Furthermore, as discussed 
in § [U a general solution to the mass discrepancy for all Cepheids is likely to require a 
combination of both mass-loss and overshoot. However, since only the mass loss can be 
observed directly, future efforts to identify large-scale circumstellar debris around a larger 
sample of Cepheids and a comparison of these results with theoretical models (including 
evolutionary models and non-linear convective hydrodynamic models) are likely to provide 
powerful new constraints on the relative contributions of mass loss versus overshoot for 
Cepheids of various masses. 

5. Conclusions 

We have used the VLA to search for Hi 21-cm line emission in the circumstellar en- 
vironment of the archetype of Cepheid variables, 5 Cephei. We have detected an extended 
(13', or ~1 pc) nebula at the position of the star. The nebula exhibits a head-tail morphol- 
ogy, with the head of the structure aligning closely with the infrared-emitting nebula and 
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bow shock previously detected by Marengo et al. (2010b), while the tail appears to be a 
turbulent structure that trails the motion of the star through the ISM. We measure an H I 
mass for the nebula of Mm ~ 0.07 M , although its total Hi mass could be ~2-3 times 
larger, depending on the fraction of emission that is hidden by the strong Galactic emission 
that contaminates a portion of our observing band. Additional material may also be present 
on large angular scales ( >15') to which the VLA is insensitive. We interpret the bulk of 
the H I nebula surrounding 5 Cephei as arising from a stellar wind with a mass-loss rate of 
M pa (1.0 ± 0.8) x 10~ 6 M Q yr~ x and an outflow velocity V Q w 35 km s _1 . By comput- 
ing evolutionary models that include mass loss across the instability strip, we show that a 
mass-loss rate of this magnitude, sustained over the preceding Cepheid lifetime of 5 Cephei, 
could be sufficient to resolve a significant fraction of the longstanding discrepancy between 
the masses of the star derived from stellar pulsation versus stellar evolutionary models. 

It is a pleasure to thank A. Pietrinferni for the computation of several evolutionary tracks 
and for many insightful discussions concerning the evolutionary properties of intermediate- 
mass stars. We are also grateful to H. Neilson for valuable discussions concerning Cepheid 
mass loss. The observations presented here were obtained through NRAO program AM998. 
LDM acknowledges support for this work from grant AST- 1009644 from the National Science 
Foundation. NRE acknowledges support from the Chandra X-ray Center, NASA contract 
NAS8-03060. 
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Fig. I. — Selected Hi channel maps from the Tapered VLA data (see Table 3). A star 
symbol indicates the position of 5 Cephei. Contour levels are (—12,— 6,— 3,3,6,12) x 1.83 mJy 
beam" 1 . The lowest contour is ~ 3c. The size of the synthesized beam is indicated in 
the lower corner of the upper left panel. Channels with Vlsr <13 km s _1 are dominated 
by large-scale Galactic emission that is poorly spatially sampled by the VLA, leading to a 
strongly mottled appearance. For this reason, only a few representative channels blueward 
of the stellar systemic velocity (Vlsr = —4.7 km s -1 ) are shown. 
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Fig. 2. — Single-dish Hi spectrum toward / =105°.l, b =0°5, illustrating the nature of the 
Galactic emission along the line-of-sight to 5 Cephei. The data were taken from Kalberla et 
al. 2005. The velocity range shown corresponds to our VLA observing band. The vertical 
bar indicates the systemic velocity of 5 Cephei. Note that the rms noise level in this spectrum 
(~ 0.09 K or f^O.57 Jy) is too high to permit detection of circumstellar emission even in the 
uncontaminated portion of the band. 
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Fig. 3. — Hi total intensity contours overlaid on Spitzer2Afim (left) and 70/im (right) images 
from Marengo et al. 2010b. The Hi contour levels are (l,1.4,2...11.3)xl8 Jy beam -1 m s -1 . 
The positions of 5 Cephei and its companion HD 213307 are indicated by star and plus 
symbols, respectively. The direction of space motion of 5 Cephei with respect to the ISM 
is indicated by a green arrow on the left panel. The insets in the lower left corners of each 
panel indicate the size of the VLA synthesized beam (96" x 88"). 
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Fig. 4. — Observed VLA H I spectrum toward the position of 5 Cephei, averaged over a single 
synthesized beam centered on the star (thin lines), with the best-fitting model Hi spectrum 
overplotted as a thick line (see § !4.2.1l for details). The observed spectrum was extracted from 
the "Robust+1" data cube (see Table 3). The portion of the observed spectrum indicated 
by a dotted line is heavily contaminated by Galactic emission and was excluded from the fit. 
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Fig. 5. — Theoretical H-R diagrams for a 5M Q star with a helium abundance V=0.273 and 
metal abundance zT=0.02. The dashed lines indicate the evolutionary path of the star. The 
solid lines depict the boundaries of the instability strip. The top panel shows a canonical 
model with no mass loss. The center panel includes mass loss starting on the main sequence 
and following a Reimers law with rj =0.4 (Eq. 2). The lower panel shows a model where 
mass loss with a rate M = 5 x M & yr -1 turns on at the red edge of the instability strip 
and continues throughout the Cepheid phase. The total time spent on the instability strip 
for each of the models is indicated in the lower right corner. No overshoot is included. See 
§ 14.51 for further details. 
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Fig. 6. — As for Figure [5J but for an initial stellar mass of 6M . 



